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We analyse the combined effect of orbital and Pauli depairing on the superconducting state, and 
apply the results to the heavy fermion CeCoIns. We find that: a) standard extrapolation based 
on the slope of H C 2(T) in the vicinity of the transition temperature to T = does not always give 
accurate values of the orbital upper critical field; b) critical value of the Maki parameter, a, that 
determines onset of the first order transition depends on the Fermi surface shape and the symmetry 
of the gap, and is a* « 3 for CeCoIns; c) the anisotropy of the thermodynamic and transport 
coefficients in the low-temperature, low-field part of the phase diagram is essentially insensitive to 
the Zeeman field and can be used to determine the nodal directions in Pauli-limited superconductors. 
The latter result confirms the finding of the d x 2_ y 2 order parameter CeCoIns. 
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Keywords: d-wave superconductors, mixed state, paramagnetic, heat conductivity, heat capacity 



I. INTRODUCTION 

Magnetic field is one of the most widely utilized and 
powerful probes of the unconventional superconducting 
state. It couples to momentum and spin of the quasipar- 
ticles, via the orbital and Zeeman mechanisms respec- 
tively, and therefore can be used to probe both the mo- 
mentum dependence of the order parameter and the spin 
structure of the Cooper pairs. 

In singlet superconductors both mechanisms are detri- 
mental to superconductivity. The orbital coupling of the 
Cooper pair motion to the vector potential of the field 
leads, in type-II superconductors, to the appearance of a 
mixed state with partial field penetration. The Cooper 
pair supercurrents and the order parameter are spatially 
modulated and form an array of Abrikosov vortices. The 
upper critical field at which the second order transition 
into the normal state occurs is given H°2 b ~ $o/27r£ 2 , 
where $o is the flux quantum, and £ is the superconduct- 
ing coherence length^. For unconventional superconduc- 
tors, below that field the number of the quasiparticlcs 
outside of the superconducting condensate, and, conse- 
quently, entropy-sensitive properties such as the specific 
heat or thermal conductivity, depend on the relative ori- 
entation of the field with respect to zeroes (nodes) or deep 
minima in the energy gap 2 - - — . The corresponding mea- 
surements have been extensively used to determine the 
symmetry of the superconducting stated—, with the 
heavy fermion CeCoIn5 and Celrln5 two examples where 
a reasonably detailed comparison of theory and experi- 
ment was carried out^ r 11 ' 15 ' 21 . Crucially, the anisotropic 
signal changes sign depending on the values of T and 
H in the superconducting phase diagram, so that cither 
minima or maxima correspond to the direction of the 
field along the nodes^— . Consequently, a detailed the- 
oretical analysis is required for the interpretation of the 



experimental data. So far, such an analysis has only been 
carried out for purely orbital coupling. 

At the same time measurements on CeCoIn5 clearly 
show that the second suppression mechanism, depair- 
ing due to Pauli spin polarization, is important in this 
compounder—. From the theory perspective, magnetic 
field aligns the spins of the unpaired electrons, while 
the singlet Cooper pairs cannot take advantage of the 
lower energy offered by a spin-polarized state. As a re- 
sult, in the absence of the orbital effects, the normal 
state is energetically favorable above the Pauli critical 
field Hp(T = 0) = A s /y/2fiB (for isotropic supercon- 
ductors, see below), where A s is the superconducting 
gap, and [ip = he/2m e c is the Bohr magneto n 29 : 30 . 
The normal state is reached via a second order transi- 



tion at T > T$ 



0.56 T c and a first order transition 



at low temperatures < T < Tp. In the absence of 
strong anisotropy of the spin-orbit coupling, the number 
of unpaired electrons at H < Hp does not depend on 
the field orientation and therefore whenever this mecha- 
nism is dominant it is not immediately obvious how effi- 
cient the field is for determination of the gap symmetry, 
and when exactly the inversion of the anisotropy pattern 
takes place. 

Under most circumstances the orbital suppression 
mechanism dominates, H°^ <C Hp. Pauli limiting is 
important in thin films and other quasi-two-dimensional 
materials, when the orbital coupling is inefficient 3 -^—. 
In strongly layered materials even at relatively low fields 
the Zeeman pairbreaking effect may complicate the ex- 
traction of the nodal contribution to the specific heat 3 -^. 
Even in relatively 3D systems, if the characteristic elec- 
tron velocity (Fermi velocity) is low and hence the coher- 
ence length is short, as it is in many heavy fermion and 
other correlated superconductors, the orbital upper crit- 
ical field is high, and can be comparable to the Pauli 
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limiting field. The question of what the vortex state 
anisotropics can tell us about the superconducting gap 
symmetry in this situation has not been explored before 
and is the subject of this paper. 

We consider a model relevant to CeCoIn 5 ) 26 ' 27 i 35 . The 
situation when the orbital and Pauli pairbreaking mech- 
anisms are comparable was considered early on by Gru- 
enberg and Gunther— and by Maki2I for an s-wavc su- 
perconductor with spherical Fermi surface (FS). CeCoIn 5 
has (i-wave gap symmetry with the main /-electron con- 
taining Fermi surface sheet open along the c-axis^Sr— . 
Consequently, as our first task, we compute the upper 
critical field as a function of temperature for different 
orientations. We extend previous calculations 4 ^— to the 
Fermi surface in the shape of a corrugated cylinder that 
gives the correct normal state resistivity anisotropy, and 
find the the temperature T*, below which the normal to 
superconducting transition becomes first-order. One of 
our main findings is that the commonly used criterion 
for determining the dominant Pauli limiting regime via 
the Maki parameter—, a = \/2H°2 b /Hp > 1 is quantita- 
tively incorrect and depends on the symmetry of the gap 
and on the shape of the Fermi surface. 

Second, we analyse the anisotropy of the specific heat 
and the thermal conductivity under rotated magnetic 
field across the T-H phase diagram including the effects 
of Pauli limiting. Since the precise details of the behav- 
ior of the anisotropy are used to infer the gap symmetry 
from experiment, this extension is critical for justifying 
d x 2_ y 2 symmetry in CeColns. We show that even mod- 
erately high Pauli effect has little influence on the the 
low temperature and low field part of the phase diagram 
whereas its consequences near the upper critical field are 
considerable. 



II. QUASICLASSICAL FORMULATION 

A. Basic equations for superconductors under 
Zeeman and orbital field. 

We use the quasiclassical formalis m 45 ! 46 at real fre- 
quencies, e, which allows to carry out calculations for 
arbitrary temperature and field and to self-consistently 
include effects of the field and impurities on the order 
parameter. The quasiclassical transport equation for the 
matrix Green's function in particle-hole and spin space, 
g, has the form 4 ^ 

[(e + ^v/(p)A(R)) r 3 - /xB • S - A(R, p) - a imp (H; e), 
g(R, p; e)] + iv f (p) ■ V R g(R, p; e) = . (1 

Here [Ei,E2\ denotes a commutator, we carried out the 
standard 4 ^ separation of the center of mass coordinate, 
R, and the momentum of the relative motion, p, so that 
the Fermi velocity depends on the position on the Fermi 
surface, v/(p). The orbital coupling is via the the elec- 
tron charge and the vector potential, A(R), and the Zee- 



man term, /iB • S, is proportional to the electron's mag- 
netic moment [i = (g/2)yUs with g-factor as a material- 
specific parameter, and the spin matrix, 

< 2 > 

with er the usual vector of Pauli matrices. The (retarded) 
Green's function in particle-hole and spin space, 

9 -{icTyif + f'cr) -g + go-* ) ' W 

satisfies normlization g 2 = —ir 2 l. Eq. ([T]) is comple- 
mented by two other equations. One is the used to de- 
termine the impurity self-energy, which we treat in the 
self-consistent t-matrix approximation 4 ^, ai mp (B,;e) = 
n imp t{R\e), with 

*(R;e) =ui+uN (g{R,p;e)}pt, (4) 

where iVo is the density of states at the Fermi level. 
In writing this equation we assumed non-magnetic 
isotropic scattering with the individual impurity poten- 
tial strength u, so that the impurity self energy does not 
have any momentum dependence 4 ^. The second equation 
is the self-consistency condition on the order parameter 
which relates it to the off-diagonal, in particle-hole space, 
component of the Green's function. Before we write it ex- 
plicitly, however, it is convenient to rewrite the function 
j in a different representation. 

Recall that we are interested in the regime H 3> H c \, 
when the vortices form an Abrikosov lattice, and are con- 
sidering strongly type-II superconductors. In this case 
the internal field B is essentially uniform and equal to 
the applied field, H. The important point is that for a 
field which orientation is the same at all spatial points 
we can choose the spin quantization axis along the field 
H = Hi, and introduce for all vector quantities the no- 
tation x ■ er = xa z . Then the 4x4 Green's function, 
Eq. ([TJ consists of two blocks, 

/ 9 + 9 / + / \ 

d R_ g-g -f + f 

9 ~ /'-/' -g + g ' W 

V -/'-/' -g-g J 

corresponding to the spin-up and spin-down components. 
The impurity self-energy 3i mp assumes the same block- 
form. 

The block structure allows us to rewrite the equa- 
tions for the components of the Green's function in a 
simple form. We explicitly introduce the equations for 
the two spin components, s = {t (+1)j4- (~ y i a 
g s = g + sg, f s = f + sf. These functions now satisfy 
independent quasiclassical equations and normalization 
conditions, g 2 — / a / a ' = — 7r 2 , with the Zeeman energy shift 
e — > e =F fiB for up and down spin respectively, see also 
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Ref. 15CJ. Now, for example, equation for the off-diagonal properly rescaled Fermi velocity component perpendicu- 



part of the Green's function takes the form 

2e 



-2i(i a -s(xB)+v f (p) V 



/.(R,p;e) 



2A 8 (R;£) l53 (R,p; £ ) 



Here /we introduced the shorthand notations e s = e — T, s 
and A s = A + A imps explicitly including the diagonal 



(S s ) and off-diagonal, Ai mp ^ components of ai mp . The 
spins mix only through the self-consistency equation for 
the singlet order parameter, 



-f-OO 

A(R ' k;£) = /^ tanh ^ 



(7) 



x ^(k,p)(/ t (R,p; £ ) + / i (R,p;e)^, 

where V^(k, p) is the pairing potential. Eq. (0 has the 
same form as the quasiclassical equation in the absence 
of Zeeman term, and therefore we can utilize the existing 
techniques to solve for each of the two spin components 
independently, enforcing the self-consistency for the or- 
der parameter at the final step. 



B. Model and method of solution 

We follow the approach we developed earlier— ~ n ' 51 
and solve the equations using a modified Brandt-Pesch- 
Tewordt (BPT) approximatio n 52 ' 53 . We assume the ex- 
istence of the Abrikosov vortex lattice and model the 
spatial dependence of the order parameter by A(R, p) = 
A(p)(R|0) = A^(p)(R|0), where y(p) is the nor- 
malized ((y(p) 2 )p = 1) basis function for the irreducible 
representation corresponding to the chosen gap symme- 
try, and normalized by a proper choice of Ck y coefficients 
the spatial vortex lattice profile, 



:r|o)=e^v- 



ikyy/SfV 



The magnetic length A 2 = hc/2eB, and 



(8) 



$oO,fcy) = *C 



A\/S7fc. 



A^/S 



/ 



(9) 



is the ground state oscillator function, and x and y are 
in the direction normal to the field.— Here we approx- 
imated the vortex lattice by only the superposition of 
the lowest oscillator wave functions, | 0). The admixture 
of the higher oscillator states is smal l 54 ' 55 and does not 
substantially affect the conclusions regarding the prop- 
erties in the vortex state d 11 ' 55 . Hence we consider only 
the lowest Landau level, as reflected in Eq. © , but with 



lar to the field. For a Fermi surface rotationally invariant 
around the c axis that we consider below^ 



(6) 



cos oh 



oil • 2n 
— sin Oh 

ox 



1/2 



(10) 



9h is the angle between the field direction and the c-axis, 



2CV 2 (p> 2 (p z )) FS , and «g ± = 2CV 2 (p>i l (p 2 )) FS , 
'II is the c-axis component of the Fermi velocity, 



o|| 

where vu 

while v±i with i = a, b is the Fermi velocity component in 
the a-b plane. For the field in the basal plane Bh = 7r/2, 
and therefore Sf = «o||/ w ox- 

The BPT approximation consists of replacing the di- 
agonal part of the Green's function with its spatial aver- 
age, and is justified over a wide range of ficld a 56 ' 57 . The 
equations for the off-diagonal components of the Green's 
function are solved by introducing the ladder operators 
as in Ref. [TolfTll and, in conjunction with the normaliza- 
tion condition, for the lowest Landau level give 



f 20FA 

fs = 7T — _ , , W 



s/T+P \v 
—i 



f 



2(e- Sf iB)A 



\vf\ 



A. s 



9s = 7T 



P 



(ii) 




w 



2(e-sfxB)A 



\vf\ 



A S A S , 



with A s (p,e) = A(p) + A imPtS (e), W{z) 
cxp(— 2 2 )erfc(— iz), and 



2f£ + »?,„ ( P) S , 



1/2 



(12) 



This closed form solution is used to enforce the self- 
consistency on the impurity self-energy and the gap 
value. The approach gives the order parameter and the 
Green's function that can be used to determine the phys- 
ical properties below. 



III. FERMI SURFACE, UPPER CRITICAL 
FIELD, AND MAKI PARAMETER 

We consider a model Fermi surface that approximates 
the main sheet of CeCoIns as seen by the magnetic os- 
cillations — ~ 41 ' 58 . It has the shape of an open cylinder 
rotationally symmetric in the a-b plane (component of 
the momentum labeled p r ) and modulated along the c- 
axis— , and is given by pj = Pr — r 2 pj cos(2ap c /r 2 pf) with 
r = a = 0.5. This choice of parameters gives a moder- 
ate anisotropy between transport coefficients in the c- 
direction and a-b plane, close to that of CeCoIns in the 
normal stated. With this definition, the typical quasi- 
particle velocity in the basal plane is Vf = pf/m, and 
that along the c-axis, Vf tC ~ avf = 0.5vf. 
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We choose the d x 2_ y 2 symmetry of the order pa- 
rameter and take a model separable pairing interac- 
tion V(k,p) = V y(k)y(p) = V y(4>W), where <f> 
labels the azimuthal angle around the Fermi surface, 
and y{<j>) = v/2cos2</>. With this choice A(R,p) = 
\/2A(R) cos 2(j). The dimensionlcss coupling constant 
NqVq determines the transition temperature of the pure 
sample, T c q, which is suppressed by impurities to T c ; we 
use this latter T c as a unit of energy. Similarly, the natu- 
ral unit for the magnetic field is the characteristic orbital 
scale, -Bo = 3>o/27r£oi wriere $o = hc/2e is the flux quan- 
tum, and £o = hvf /2itT c is the coherence length. 

We measure the strength of the Zeeman term via a di- 
mensionless parameter Z = /.iBq/2ttT c ^- so that the Zee- 
man splitting of the energy levels in dimensionless units 
is \iB I (2irT c ) = Z B/Bq. Since we expect the orbital 
critical field H°^ b to be of order Bq, and the Pauli limit 



(a) tf™*=1.53B (llantinode) 
H_H " 1 



(b) H f= 0.57 B (II c-axis) 



to correspond to fj,H P ~ A ~ T c , we find H° c2 ~ b /H P ~ Z. 

We first compute the upper critical field by solving the 
linearized, with respect to the order parameter, quasiclas- 
sical equations. Under the combined effect of the orbital 
and Zeeman field the transition for strong enough Pauli 
term becomes first order at low temperatures, T < T* . 
We determine T* in the clean limit by evaluating where 
the A 4 -term coefficient in the free energy expansion be- 
comes negative. The general free energy expansion can be 
obtained from Eqs. but it is complicated in a dirty 

d-wave superconductor—, and since the precise location 
of T* is not important for this work we will not look for 
it here. At lower temperatures, the transition line can 
only be determined from the full free energy functional, 
which needs to also account for the possible existence of 
the additional modulations in the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase^r— . While there are many 
indications that in CeCoIns a new phase exists in the 
high-field, low-temperature rang e 35 ' 62 " — the experiments 
on the vortex state anisotropy are generally carried out 
away from that range. Consequently, we do not consider 
the FFLO-like modulation here. 

Fig. [T] shows the upper critical field, H C 2(T) for a 
pure system for two different field orientations (c-axis 
and the in-plane along the gap maximum) and for dif- 
ferent strength of the Zeeman splitting. Our results 
for Z = resemble the classical results by Helfand 
and Wcrthame r 67 ' 68 (HW) with differences in profiles at- 
tributed to anisotropic gap and non-spherical FS. The 
important observation is that the shape of the tempera- 
ture dependence of the upper critical field changes with 
increasing Z, and that the linear region of H C 2(T) ~ 
1 — T/T c near T c rapidly shrinks, eventually leading to 
\J\ — T/T c dependence characteristic of the Pauli 



H, 



limited field. Hence, as is shown in Fig. [T] a frequently 
used experimental estimate of the H°2 b (T = 0) based 
on the slope of the measured upper critical field near 
T c is not reliable for paramagnetically limited supercon- 
ductors: that estimate changes on approaching the Pauli 
limit. First, the slopes dH C 2/dT of the purely orbital and 
paramagnetically limited cases, while equal asymptoti- 




FIG. 1: (Color online) Illustration of the possible pitfalls 
in the experimental determination of orbital limit on upper 
critical field, H°2 b , in strongly paramagnetic superconduc- 
tors, (a) in-plane field, (b) field normal to the planes. Lower 
(blue) solid lines denote H C 2 for the paramagnetically limit- 
ing case Z — 0.8, and the dashed lines mark regions of the 
first-order transition. "Turning off" Zeeman effect (Z = 0) 
results in true orbital H C 2, the upper (black) solid lines. Dot- 
ted lines: slope dH C 2/dT at T c determined from the points 
T/Tc — (0.9, 1). Thick dot-dashed lines: upper critical field 
obtained using the Helfand- Werthamer (HW) result for the 
given slope dH C 2/dT at T c . Thin dot-dashed lines: HW pro- 
file of H C 2{T) for the paramagnetically limited case Z = 0.8, 
resulting in incorrect extracted orbital limit, H%2. Vertical 
arrows indicate potential discrepancy between the approxi- 
mately determined and realistic orbital upper critical field, 
see text for details. 



cally as T — > T c , are different when determined within a 
reasonable experimental window as shown by dotted lines 
based on the values of H c2 at T/T c = (0.9, 1). We find 
that even for moderate paramagnetic coupling a reason- 
able estimate of the slope may be obtained only within a 
window of 1-3% near T c . Second, the H c2 {T) profile for 
our chosen Fermi surface and the d-wave order parame- 
ter even for the purely orbital case is different from that 
found using the HW— result for s-wave superconductor 
with a spherical FS with the same dH&jdT (T c ). Con- 
sequently, for Z ~ 1 the estimate of H°2 b using the HW 
profile and an approximate slope dH C 2 j dT (T c ) based on 
points more than a few percent away from T c (thin dot- 
dashed lines) may significantly underestimate the mag- 
nitude of H°2 b , as indicated by the vertical arrows in 
Fig.ffl 

The change in the shape of the H C 2(T) curve with in- 
creasing paramagnetic contribution is made even more 
explicit in Fig. [Jl where we also include the upper criti- 
cal field for the in-plane field along the nodal direction. 
As can be expected in a magnetically isotropic system, 
increased Pauli limiting reduces the anisotropy of the up- 
per critical field caused by both the in-plane anisotropy 
of the order parameter (node vs. anti-node) and the the 
anisotropy of the Fermi surface (c-axis vs. a-b plane). 

Finally, in the inset of Fig. [5] we show the onset tern- 
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FIG. 2: (Color online) Upper critical field of a pure d x 2_ y 2 su- 
perconductor with quasi-cylindrical FS for three orientations 
of the field. The critical value of the Maki parameter is about 
a* ~ 3 for all directions (see text). The dotted lines indicate 
first-order transition below T* . Inset in the right-hand panel: 
T* as a function of the relative strength of the Pauli limiting. 
For purely Pauli-limited case Tp ~ 0.56T C . 



perature for the first order transition, T* as a function of 
the Pauli limiting parameter Z. It is instructive to recast 
this analysis in terms of the so-called Maki parameter, 
a = \/2H°2 b I Hp. It is conventionally assumed that the 
critical value of the Maki parameter that defines a strong 
paramagnetic limit is a > a* = 1, as obtained for the on- 
set of the first order transition in s-wave superconductors 
with a spherical Fermi surface 3 -^; this value is commonly 
used in the analysis of anisotropic strongly correlated ma- 
terials as well. In reality the critical value a* depends 
on details of the Fermi surface and the superconducting 
state. For our model H°2 h ~ 1.4 in the a-b plane, and 
\iH P ps l.HA d o/\/2. Here A d0 = 0.241(2ttT c ) is the am- 
plitude of the order parameter at T = H = 0— and the 
increase of Hp by factor 1.11 compared with the simi- 
lar s-wavc expression is caused by the gain in magnetic 
energy of the d-wave superconducting state by magneti- 
zation of the nodal quasiparticlesSi. The critical value of 
the Maki parameter where first order transition appears 
is therefore a* = V2H% b /H P ps 7.0(H^ b / B )Z* , and, 
since in-plane Z* ~ 0.3 we have ot* d _ wave ps 3.0, signif- 
icantly exceeding the value of unity expected for dirty 
isotropic systems. 

We apply our results to CeCoIn 5 to analyse the be- 
havior of the upper critical field along a and c axes. The 
fit for the [100] direction is shown in Fig. [3] The best 
fit of the second order transition and the location of T* 
in the experimental range is given by Z a ~ 0.5. For 
H\\c the best fit is given by Z c ~ 1.15. From these two 
fits we obtain the value for Bq which is approximately 
30 Tesla for both fields orientations, which is an indica- 
tion that our choice of the FS parameters was reasonably 
good to describe this system. From this value we find 
the Fermi velocity Vf ~ 0.6 • 10 6 cm/s, which is sensible 
for a heavy fermion system and agrees with interpreta- 



FIG. 3: (Color online) A fit of the experimental H c2 // [100] in 
CeCoIns. The data are taken from measurements of magneti- 
zation by T. Tayama et aii— and specific heat by A. Bianchi 
et a/.—. We find that Z — 0.5 gives the best fit of second 
order transition line and the onset of the first order transition 
T*. 



tion of the experimental data^, and disagrees with the 
value ~ 10 8 cm/s obtained by authors of Refl44l. We 
then extract the values of the effective electron moment, 
m/mb = Z{2itk,BT c ) / (ijlbBo) and obtain u a ,b/p<B ~ 0.35 
and Hc/l^B ~ 0.7, in agreement with Ref|44j where the 
varying parameters were Fermi velocity and the g-factor. 

We do note that in the current model the orbital pair- 
breaking is stronger for the nodal direction, Fig. HJa), 
hence the Pauli effect is relatively less important, and the 
range of first order transition is smaller for that orienta- 
tion, (Figrjjc)). However, this particular aspect depends 
sensitively on the in-plane shape of the Fermi surface. We 
took the Fermi surface to be rotationally symmetric, and 
the question of how this conclusion is affected by a more 
realistic Fermi surface shape is left for future studies. 



IV. TRANSPORT AND SPECIFIC HEAT 
ANISOTROPY 

We are now in the position to analyse how the Pauli 
limiting affects the measured anisotropies of the spe- 
cific heat, C(4>o), and longitudinal thermal conductivity, 
K xx ((f>o), when magnetic field is rotated in the basal plane 
with respect to the crystalline axes. The direction of the 
field, (f>o, is measured from the gap's maximum along the 
a-axis. As is well known, both quantities show oscilla- 
tions, with either minima or maxima when the field is 
aligned with the nodal directions^ -11 ' 22 ' 23 . The loca- 
tion of the inversion line in the T-H plane that separates 
the regions with minima and maxima indicating the gap 
nodes is of exceptional experimental relevance as it af- 
fects the conclusions about the shape of the gap in a 
given compound. 

The results are most clearly presented in the form of a 
phase diagram in the T-H variables that indicates each 



6 



H/B tj =0.1 {circled 



WB =0.25 (squares) 



03 
5 



Z=0.0 



.1.1.1.1 




- (el) 


7/T 

0.6 
0.5 
0.4 
0.3 


- (c2) - 








0.15 
0.1 






0.05 






0.2 0.4 0.6 0J 
T/T 



1 0.2 0.4 0.6 O.t 
T/T 



J 0.2 0.4 0.6 0.8 
T/T 



FIG. 4: (Color online) The anisotropy of heat capacity phase 
diagram for three Zeeman couplings Z = 0.0 (a), Z = 0.2 (b) 
and Z = 0.4 (c). Shaded areas correspond to minimum of the 
heat capacity for H|| (node). Panels (cl) and (c2) demonstrate 
how anisotropy of heat capacity, C(<j>o), varies with tempera- 
ture and magnetic field for Z — 0.4. The curves for different 
T/T c are shifted vertically for clarity. 



of the regions. In Fig. [4] we show the evolution of such 
phase diagram for the specific heat anisotropy with in- 
creasing Zeeman coupling. We compute the density of 
states, N(e,T) = — (N Q /2ir)(3mg(p,e))p, where inclu- 
sion of T signifies that the self-consistently determined 
gap is temperature-dependent, obtain the entropy, 



/oo 
dojN(e, T) 
-OO 



/(e) In /(e) 



+ (l-/( £ ))ln(l -/(e)) 



,(13) 



and numerically differentiate it to find C/T. The result 
for purely orbital coupling to ma gnet ic field, Z = 0.0 in 
Figffla,) is in agreement with Ref. [HJ In the two shaded 
regions (near H C 2 and in low-T low-H corner) the heat 
capacity attains its minimum when the field is along the 
nodal directions of the gap. Clearly, since the anisotropy 
in the heat capacity is field-induced, the amplitude of the 
anisotropy is very small at low fields for all temperatures, 
and extremely small near T c . Hence the challenge is to go 
to sufficiently high fields to see the signal while knowing 
whether the maximum or a minimum of the anisotropic 
signal corresponds to the nodal directions. 

As the relative strength of the spin coupling increases, 
Fig. SJb) and (c), we note that the "nodal minimum" 
region near iJ C 2 grows at the expense of the intermediate 
field unshaded region. This is concomitant with the dis- 
appearance of the anisotropy of the upper critical fields 
for nodal and antinodal directions. An obvious conjec- 
ture is that for Z = at intermediate temperatures and 
moderately high fields it is the anisotropy of the upper 
critical field that controls the anisotropy of the specific 

, at a fixed external field 
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< H; 
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FIG. 5: (Color online) Anisotropy of the specific heat (left) 
and heat conductivity (right) for Z = 0.0 and Z = 0.4 at 
T/T c — 0.3 and various fields. The field values are shown on 
the left in terms of Z-dependent H C 2- The curves were shifted 
vertically for clarity. The absolute values of C and k and their 
anisotropics presented in Figures below. The range of fields 
in H/Bq for convenience indicated on the right. 



with the heat capacity, is higher for the field along the 
nodal direction. As the Zeeman coupling is increased and 
the critical fields along the nodal and the antinodal direc- 
tions become nearly equal, this unshaded region shrinks. 
In Fig. @Jc) we tuned Z > Z* , and see that the criti- 
cal fields for the two in-plane directions are almost in- 
distinguishable, and that there is a region of first order 
transition at low temperatures T < T* , indicated by the 
arrows. 

Importantly for our purposes, the low-T, low-H shaded 
region (minima for the field along the nodes) remains 
mostly unchanged with increasing Z, terminating at 
T/T c - 0.1 and H/H^ ~ 0.4. This region is still dom- 
inated by the nodal quasiparticles and "semiclassical" 
physics, A(0, H) - y/W/H^P^. compared to a lin- 
ear contributions due to Zeeman shifts (for these Z's) 
and vortex cores. 

The insets Figs. BJcI) and@Jc2) provide image of the 
heat capacity anisotropy for Z = 0.4 coupling for two dif- 
ferent fields at locations indicated by circles and squares 
correspondingly. The heat capacity expression that we 
use to plot the angular dependence, 



C(T,H) 



de e 2 N t (e, T, H) + N±(e, T, H) 



4T 2 



cosh 2 (e/2T) 



(14) 



H / H node > jjjjjantinode^ and the d ensity f States, along 



is, strictly speaking, valid only at low temperature 
when the order parameter is essentially temperature- 
independent. The difference between this expression and 
the exact result obtained from the entropy differentia- 
tion is small far from the superconducting transition, see 
Refs. HH 

We also compute the anisotropy of the thermal con- 
ductivity along the crystalline x direction for unitarity 
scattering (phase shift it/ 2) and the normal state scat- 
tering rate T/2itT c = 0.007. The thermal conductivity 
is calculated on equal footing with the density of states 
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FIG. 6: (Color online) Co</> and coefficients in heat ca- 
pacity expansion at two different temperatures, T = 0.05T C 
and T = 0.3T C . Positive four-fold coefficient means that the 
specific heat has minimum for the field along a node in the 
gap, while negative C44, corresponds to the maximum for the 
field along a node. 



and is given by^ 



11.51 



k xx {T,H) f 

2T T 2 



de e 2 , _ 2 e 
^cosh — 



(15) 



x (vj iX N(T,H;p,s)T H (T,H;p,e] 



with the effective scattering rate 



2t h 



= -3m E 



R 



2A 3m[g R W{2eA/\vj\ 



/ 



Jm g 1 



|Ao^| 2 



(16) 

Here the label R denotes a retarded function, and we use 
the angle-resolved density of states. 

Fig. [3] shows characteristic profiles of the heat capac- 
ity and heat conductivity when the field is rotated in the 
basal plane. Note that in the left panel we labeled the 
fields according to the ratio H/H C 2. the upper critical 
field changes between Z = and Z = 0.4, so that the 
effective field range is different. As discussed above, the 
major difference between two values of the Zeeman split- 
ting is that for higher Z there exists a region of shallow 
minimum at the node (45° in our case) for high fields. 

The second panel of Fig. [5] shows the behavior of the 
thermal conductivity for the same fields, now labeled 
in units of Bo. Qualitatively, the peak in the angle- 
dependent k for Z = at 45° disappears at Z = 0.4, 
making the dependence of the thermal conductivity more 
twofold. 

To quantify these trends we follow the approaches 
taken in experiment and expand both quantities in har- 
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FIG. 7: (Color online) The anisotropic component of the 
specific heat at temperatures T — 0.05T C (left panel) and 
T = 0.3T C (right panel) for Z = and Z = 0.4. 



monies of the angle between the field and the heat cur- 
rent, 



c(<fe,r,g)/r 

(C/T) N 

{Kxx/T)n 



= C 0(f> + C 44> cos4:(l) Ql (17) 

= K O -(- K 2c j, COS 200 + ^40 cos 40o , 



and present the evolution of different coefficients with 
field and temperature. 

First, focus on the behavior of the specific heat as 
shown in Fig. [6] At low fields for Z = the isotropic 
part, Co(f>, shown in panels (a) and (c), exhibits the ap- 
proximate v^ff behavior expected of a nodal supercon- 
ductor— With increased Zeeman contribution this 
component of the specific heat becomes more linear in 
g e i c j3 i 75_ Thj s j Sj f coursej simply because the spin-split 
quasiparticle spectra produce a density of states at zero 
energy, iV(0, H) ~ (fj,H) /Aq, with the prefactor that, for 
high n, can be large enough to dominate the sublincar 
Volovik term already at low fields. At low temperatures 
Co0 acquires a pos itive curvature at high fields, in agree- 
ment with Rcf. [6fJ. 

The qualitative behavior of the anisotropic term,, 64^, 
is similar for the two values of Z at low temperature, 
panel (b), but is distinctly smaller at higher tempera- 
ture, shown in panel (d), and even changes sign near 
H C 2- The comparison becomes even more clear if the 
amplitude C40 is plotted against the upper critical field 
for given values of Zeeman splitting, as in Fig. [71 At low 
temperature there is essentially no difference between the 
purely orbital case and that of a moderately strong Zee- 
man coupling, and hence in this regime the anisotropy 
depends only on the shape of the gap and the underly- 
ing Fermi surface. In contrast, at T/T c — 0.3 for the 
case of weak Zeeman splitting the coefficient does not 
change sign, as seen already from Fig. QJa). In contrast, 
for Z = 0.4 as shown in Fig. Stc), there is a change of 
sign in the anisotropy of the specific heat, due to the 
"isotropization" of the upper critical field, as discussed 
above. 

The same trend is seen in the thermal transport 
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FIG. 8: (Color online) Field dependence of the coefficients Ko 
and K2 in the expansion of the thermal conductivity, Eq. (|I8[) 
at T = 0.05T C (top two panels) and T = 0.3T C (bottom pan- 
els) for Z = and Z = 0.4. 



anisotropy. The behavior of the average thermal con- 
ductivity, kq and the twofold term responsible for the 
difference between heat transport parallel and perpen- 
dicular to the vortices, is shown in Fig. [8l At low fields 
the scattering of the quasiparticles on the vortices is de- 
termined by the vortex concentration, n ~ ijy$o, an d 
therefore, at low temperature, only weakly depends on 



i 1 i 1 i ■ i 1 
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FIG. 10: (Color online) The phase diagram for the anisotropy 
of the specific heat under rotated magnetic field for different 
strength of Pauli pairbreaking term. Shaded areas correspond 
to minimum of the heat capacity for H|| (node). Notice that 
the shaded region at low T and H is identical in all three 
panels. 



the Zeeman field. As the temperature increases, however, 
the differences between the two cases become much more 
pronounced. In particular, it is worth noting that the 
two-fold symmetry k 2 > (K xx (j h \\H) > K xx (j h _L H)) 
is much more pronounced in the Pauli-limitcd case. 

The fourfold "nodal" term, Fig. [SJ shows the behav- 
ior broadly similar to that of the anisotropic component 
of the specific heat. The low-temperature behavior as a 
function of the reduced field, H/H C 2 is similar for the 
cases of weak and strong Pauli limiting, while the be- 
havior at moderate temperatures is very different start- 
ing at moderate fields. Once again, the coefficients for 
Z = and Z = 0.4 have the opposite signs starting at 
H - 0.5H c2 . 




FIG. 9: (Color online) The anisotropic fourfold term of the 
thermal conductivity in Eq. (|18p is shown for T — 0.05T C in 
panels (a) and (c), and for T = 0.3T C in panels (b) and (d) as 
function of the field, panels (a) and (b) and the reduced field, 
H/H C 2, panels (c) and (d). 



V. CONCLUSIONS 

In conclusion, we calculated in a d-wave superconduc- 
tor with quasi-cylindrical FS the H C 2 and its evolution for 
transition between orbital and Pauli limits. We find that 
in this system the critical value of the Maki parameter 
for paramagnetic limit criterion is a* = 3. That values 
depends on both the symmetry of the order parameter 
and the shape of the Fermi surface, and hence there is no 
universal threshold value for that parameter that ensures 
the first order transition in a given system. We also find 
that the linear extrapolation of the upper critical field in 
Pauli-limitcd superconductors from the vicinity of T c to 
low temperatures does not accurately predict the value 
of the orbital critical field. 

We showed that moderately large Zeeman contribu- 
tion does not alter the typical behavior of C— and 
k — anisotropy at low temperatures and fields. This is 
summarized in the phase diagram in Fig. 1101 The re- 
gion at low T and H where the minima in the fourfold 
terms in the heat capacity and thermal conductivity oc- 



9 



cur for the field along the nodes is essentially insensitive 
to the strength of the Zeeman coupling. The differences, 
of course, occur at higher temperatures and fields. 

This finding proves that, even in systems with strong 
paramagnetic contribution where H C 2 anisotropy in the 
a-b plane is largely absent and therefore cannot be used 
to infer the nodal structure^, the anisotropy of heat con- 
ductivity and heat capacity at moderate to low H and 
T still can help determine the location of the nodes of 
the order parameter on the Fermi surface. The previous 
analysis^ of the inversion of the anisotropy in CeCoIns 
remains valid when the Zeeman term is accounted for, 
and unequivocally indicates the d x 2_ y 2 shape of the or- 



der parameter. This also provides further support for the 
interpretation of Ref. HE] that confirmed the inversion of 
the specific heat oscillations. 
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